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This study discusses about the optimization of spray characteristics of biodiesel fuel in a spray chamber. Three 


factors namely, fuel injection pressure, fuel temperature and fuel blends were chosen as the influencing factor for 


the set objective. Four levels were chosen in each factor and spray tip penetration, spray cone angle and Sauter 


mean diameter (SMD) were taken as the response variables. Experiments were designed by employing design of 


experiments method and Taguchi full factorial array was used to conduct the tests with different levels of the 


chosen factors. Multi Response Signal to Noise ratio (MRSN) was calculated for the response variables and the 


optimum combination level of factors was obtained simultaneously using Taguchi’s parametric design. Confirma- 


tion experiments were conducted for the obtained optimum combination level of factors and the results were 


compared with the normal operating conditions and significant improvement were observed in the response va- 


riables. 
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Introduction 


Taguchi is the developer of the Taguchi method. It is a 
powerful tool for the design of high quality systems. It 
provides a simple, efficient and systematic approach to 
optimize performance, quality and cost. Taguchi design 
of experiments is most extensively used to determine the 
parameter values or setting required to achieve the de- 
sired function. Taguchi method can also be used to inves- 
tigate the effects of interactions between the various fac- 
tors which can be easily missed when using conventional 
methods. In recent years, the rapid growth of interest in 
the Taguchi method has led to numerous applications in a 
world wide range of industries and nations. Some of the 
literatures covering Taguchi approach are reviewed and 
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presented for knowing its application in engine and man- 
ufacturing related work. 

The simultaneous optimization study on diesel engine 
design and operating parameters for low exhaust emis- 
sions using Taguchi method. A single cylinder, research 
diesel engine equipped with a high pressure, cam driven, 
electronic unit injector was used for optimization expe- 
riment. The measurements obtained after optimum para- 
meter settings showed that the particulates and NO, 
emissions were significantly reduced. The Taguchi me- 
thod was used to be a useful technique for the simulta- 
neous optimization of engine parameters and exhaust 
emissions[1]. Taguchi method for optimization of critical 
peening parameters and this method is used to formulate 
the experimental layout, to establish the order of predo- 


minance among the identified critical parameters and 
predict the optimal setting for each of the process para- 
meters. The experimental results obtained were con- 
firmed the adequacy and effectiveness of Taguchi ap- 
proach[2]. To optimize the parameters such as rubber 
seed biodiesel blend, injection pressure and applied load 
of the single cylinder diesel engine with respect to spe- 
cific fuel consumption and brake thermal efficiency 
through the Taguchi and ANOVA technique method was 
used. The optimum values of the parameters were found 
using Signal-Noise ratio. Analysis of variance (ANOVA) 
was used to investigate the influence of parameters on 
the response [3]. To optimize three parameters such as 
blend proportion, injection timing and injection pressure, 
a simultaneous optimization method called Taguchi was 
used in their work. This method requires fewer numbers 
of trails for fixing optimum levels. In addition, An 
ANOVA was also performed for the parameters to eva- 
luate its percentage contribution over the desired out- 
put[4]. 

The spray characteristics of fuel play a vital role in the 
combustion process, performance and emission characte- 
ristics of diesel engine. Generally, small droplet size and 
large spread area of the fuel blends spray are considered 
as important factors to improve the atomization ability 
and to enhance the combustion efficiency. Variation of 
fuel injection pressure, fuel temperature, type of fuel 
blends has a considerable effect on spray parameters of 
fuel. Hence fuel injection pressure, fuel temperature and 
fuel blends are selected as factors for the present investi- 
gation. The higher viscosity causes mixture to burn lean 
in the engine as fuel moves slowly through fuel filters 
and fuel lines. 

In spray study, the inertia force and air drag force are 
more important factors compared to the viscous force and 
surface tension force [5].Droplet size and its distribution 
follows the vibration and breakup process. The effects of 
spray flow rate, spray height and inlet temperature on 
spray cooling were investigated and the corresponding 
droplet axial velocity and SMD were correlated with 
mean absolute error of 15% [6]. As the ratio of biodiesel 
in the blends increased, spray tip penetration increased, 
but the spray cone angle decreased [7]. A reduced fuel 
viscosity leads to a slim spray that is characterized by a 
reduced spray width and an increased spray penetration 
[8].Spray characteristics of the fuel mainly depend on 
fuel injection pressure, fuel density, fuel viscosity, am- 
bient pressure and temperature [9].The penetrating speed 
during the initial stage is primarily controlled by the 
competition between the inertia and surface tension [10]. 
Atomization and mixing of sprays are key parameters to 
successfully describe and predict combustion in di- 
rect-injection engines [11]. Spray tip penetration (spray 
length) is the distance between the injector tip and long- 


est spot in the spray image and is obtained using graphi- 
cal tool of pro analyst software. Spray cone angle is the 
angle formed by two straight lines drawn from the tip of 
the injector to the outer periphery of the spray at a dis- 
tance of one third of the spray length (s) downstream of 
injector tip and is obtained by drawing appropriate lines 
in the captured images using graphical tool. 


Experimental setup and procedures 


Fig.l shows the experimental setup for investigating 
the macroscopic spray characteristics of preheated bio- 
diesel. It includes a spray chamber, preheating setup, fuel 
injection system, high speed video camera and a data 
acquisition system. Fuel injector is mounted on the top of 
the spray chamber. The injector is a simple mechanical 
fuel injector with delivery opening pressure varying from 
180-240 bar. Fuel injector was connected to the fuel 
pump via a high pressure line and fuel pump is operated 
using an electric motor (0.25HP) for the experiments. 
The technique used for visualizing the liquid phase pene- 
tration is shadowgraphy technique with halogen lamp. 
The region of the spray has been illuminated by a halo- 
gen lamp and light scattered by liquid droplets has been 
captured by high speed video camera at the frame rate of 
1250 fps with resolution of (800x650) using Fastmotion 
software and further analyzed by ProAnalyst software. 
The ProAnalyst software was used to calculate the spray 
tip penetration(S) and spray cone angle(O) for each 
frame. The spray characteristics of tested fuel blends for 
various injection pressures and fuel temperature were 
measured by applying image processing technique, by 
frame analysis for some random selected pictures. In 
order to ensure the repeatability of the experimental re- 
sults each spray test was repeated four times and the val- 
ues given in the study are the average of four measure- 
ments. Taguchi based design of experiments (DOE) me- 
thod was employed to design the experiments for spray 
characteristics of biodiesel fuel. 


Fuel tank 


Spray chamber 


Fig. 1 Schematic diagram of experimental setup with fuel 
preheating arrangement 


Selection of levels of factors 


Factors chosen for the present investigations are not 
discrete can be measured on a scale. To find the effects of 
fuel injection pressure, fuel temperature and fuel blends 
on spray tip penetration, spray cone angle and SMD, 
their levels haven to be chosen from a minimum value to 
maximum value. For the present work, minimum and 
maximum value for the factor chosen based on the earlier 
research work conducted with those factors individually. 
Any value in between the minimum and maximum can 
be assumed to increase the levels of factors. For the 
present work, four levels were chosen for each factor to 
critically examine the effects of selected factors on the 
chosen objective as shown in Table 1. 


Table 1 Selection of factors 


Factor Factors influencing Levels of factors 
No the objective 1 2 3 4 
1 Fuel blends JB20 KB20 JOME KOME 
Injection 
2 180 200 220 240 
pressure, bar 
Fuel 

50 60 70 80 


temperature, °C 


Full factorial experimental design 


In full factorial experiment for three factors with four 
levels, the number of experiments to be conducted will 
be 4° = 64. Selection of full factorial experimental design 
matrix is shown in Table 2. In the present investigation, 
degree of freedom for the factors with four levels is 27 (3 
x 3 x 3) and degrees of freedom for Le4 are 63 (number 
of trails-1). 


Analysis of data 


Three variables (spray tip penetration, spray cone an- 
gle and SMD) were chosen as the responses of the prob- 
lem. The responses obtained for each trail at different 
operating conditions were analyzed to get a result for the 
formulated problem. In the analysis, average values of 
the responses measured at different operating conditions 
were considered as the responses for the trail. To optim- 
ize the combination of the level of the factors for the 
formulated problem, Multi Response Signal to Noise 
ratio (MRSN) was calculated. 

The procedure employed in the optimization process is 
described below 


Loss Function 


Loss function is used to calculate the deviation be- 
tween the experimental value and desired value. For each 
response variable, the corresponding loss function can be 


expressed as given below[12]. As per the Taguchi’s cate- 
gorization of response variables, smaller the better prin- 
ciple is considered to minimize the spray tip penetration 
and SMD. For spray cone angle, larger the better prin- 
ciple is considered to maximize it. 

For larger the better (Spray cone angle) 


I~ 1 
Lj = 5 (1) 
n 7 
kal Y Ük 
For smaller the better (spray tip penetration and SMD) 
1 n j 
L;==> Y n (2) 
k=1 


where n is the number of repeated experiments , Lj is 
the loss function of the im response variable in the ji, ex- 
periment at the ka test. 


Normalizing the Loss Function 


Since the measured units of the response variables 
were dissimilar, the loss function was normalized in the 
range between zero and one. Normalization of loss func- 
tion was done as follows[12]: 

For smaller the better (Spray tip penetration and SMD) 


S. = min L; (3) 
ij La 
For the larger the better (Spray cone angle) 
ç a Es (4) 
id max Lj 


Where S; is the normalized loss function for the re- 
sponse variable in jm experiment, L; is the loss function 
for the im response variable in the jm experiment and L; is 
the average loss function for the i,, response variable. 


Assigning weighting factor 

In multi response optimization, the relative signific- 
ance of each response variable on the set objective with 
respect to others will be fixed by assigning proper 
weighting factor for each of the normalized quality loss 
function by including the weighting factor, the total loss 
function (TL;) can be expressed as 


TL; = Ns) (5) 


Where w; is the weighting factor for the im response 
variable and m is the number of response. Weighting 
factors for the response variables are to be determined 
based on the priorities among the various responses. If 
equal importance is given to all the response variables, 
the weighting factors will have equal value such that the 
sum of the weighting factors is always unity. In an opti- 
mization process with three response variables, for the 
combination 0.4, 0.3 and 0.3, the importance on first re- 
sponse variable is more when compared to the other two. 


In this way different combination as per the chosen ob- 
jective can be taken to get optimum combination level of 
the influencing factors. The most influencing factor in 
achieving the objective for each permutation of the 
weighting factor was analyzed through Analysis of Va- 
riance (ANOVA). 

The main objective of the work was to reduce the 
spray tip penetration with minimum SMD and maximum 
spray cone angle. Hence higher weighting was assigned 
to spray tip penetration, when compared to the other two. 
Initially 0.4(w,), 0.3(w2) and 0.3(w3) were assigned as 
weighting factors for the response variable spray tip pe- 
netration, spray cone angle and SMD respectively. 

Further it was varied to study the effect of weighting 
factor on the set objective. 


MRSN 


In multi Response optimization of Taguchi loss func- 
tion, Multi Response Signal to Noise ratio (MRSN) has 
to be maximized by using the formula given below[12] 

MRSN =-—10log(TL ;) (6) 


Optimal level of combination for the obtained MRSN 
ratio with the assigned weighting factor was determined 
by following Taguchi parametric design. Variance of the 
MRSN ratio was analyzed through Analysis of Variance 
(ANOVA) and the level of importance of each factor on 
the response variables for the assigned weighting factor 
was identified from the ANOVA table. This procedure was 
repeated for different combinations of weighting factors 
to predict the effect of weighting factor on the set objective. 


Analysis of variance (ANOVA) 


ANOVA is a statistical method used to interpret expe- 
rimental data and make required decisions and it estab- 


lishes the comparative significance of factors in terms of 
their percentage contribution to the response. Since three 
factors are involved in the present investigation, it is ne- 
cessary to evaluate the significance and percentage con- 
tribution of each factor on the reduction of spray tip pe- 
netration with less sacrifice on spray cone angle and 
SMD. This analysis is performed on Signal to Noise ratio 
to find the contribution of the factors using Minitab-14 
software. MSF is equal to SS; divided by the number of 
degree of freedom (DF) associated with the factors. 
The F-ratio provides a statistical value that can be com- 
pared to a probability distribution table for a given con- 
fidence level to identify the significant effect of each 
influencing factor on the responses. 


Result and discussion 


Confirmation in experiment 


Optimum combination of factor levels obtained 
through MRSN ratio and Taguchi parametric design was 
confirmed through the experiment. 


MRSN ratio 


Table 2 shows the MRSN ratio for the experiments 
conducted for the weighting factors w,= 0.4, w. = 0.3, 
and w; = 0.3. From the Table 2, the combination which 
has the maximum MRSN ratio will be taken as the best 
combination among 64 in achieving the objective. It can 
be observed that experiment number (4, 4, 1) is the best 
combination among the 64 experiments and Fig. 2 shows 
the main effect plot for MRSN ratio. ANOVA was em- 
ployed to analyze the MRSN ratio obtained with different 
combinations of weighting factors. 


Table 2 Taguchi full factorial design and MRSN ratio for w, =0.4, w,=0.3, w3=0.3 


Exp full factorial Loss function (Lij) a. ae Total loss MRSN 
. Normalization (Sj) Weighting Factor : | 

No design S 0 SMD function (TL) ratio 
1 1 1 2498 0.004 113.84 0.92 0.76 0.69 0.37 0.22 0.20 0.784 0.9242 
2 1 2 2410.8 0.004 111.5 0.958 0.71 0.71 0.38 0.21 0.213 0.797 0.90209 
3 1 3 2400 0.0041 108.9 0.963 0.67 0.72 0.38 0.20 0.218 0.906 0.94297 
4 1 4 2361.9 0.0040 107.6 0.978 0.65 0.73 0.39 0.19 0.220 0.892 0.92734 
5 2 1 2554.2 0.006 102.4 0.905 1 0.77 036 0.3 0.232 0.847 0.48644 
6 2 2 2475.6 0.0055 100.3 0.933 0.89 0.78 0.37 0.26 0.236 0.833 0.55987 
7 2 3 2376.5 0.0051 98.01 0.972 0.82 0.80 0.38 0.24 0.242 0.882 0.55638 
8 2 4 2311.6 0.0044 96.7 1 0.72 0.81 0.4 0.21 0.245 0.871 0.64244 
9 3 1 3064.7 0.0051 93.04 0.754 0.82 0.85 0.30 0.24 0.255 0.759 0.9402 
10 1 3 2 2813.2 0.0044 91.18 0.821 0.72 086 0.32 0.21 0.260 0.777 0.93834 
11 1 3 3 2764.6 0.0040 89.05 0.836 0.65 0.88 0.33 0.19 0.266 0.824 0.97734 
12 1 3 4 2621.4 0.0037 86.5 0.881 0.60 0.91 0.35 0.18 0.274 0.829 0.91821 
13 1 4 1 3137.1 0.0044 85.2 0.736 0.71 0.92 0.29 0.21 0.27 0.752 1.02741 
14 1 4 2 3047.0 0.0041 83.53 0.758 0.66 0.94 0.30 0.19 0.284 0.762 1.04510 
15 1 4 3 2916 0.0035 81.28 0.792 057 0.97 0.31 0.17 0.292 0.844 1.06578 


Continued Table 2 


Exp full factorial Loss function (Lij) o a Total loss MRSN 
. Normalization (Sj) Weighting Factor . á 
No design S 0 SMD function (TL;) ratio 
16 1 4 4 2787.8 0.0034 79.21 0.829 0.55 1 0.33 0.16 0.3 0.723 0.97679 
17 2 1 1 2809 0.0047 132.5 0.822 0.75 0.59 0.32 0.22 0.179 0.721 1.32925 
18 2 1 2 2652.2 0.004 122.9 0.871 0.70 0.64 034 0.21 0.193 0.741 1.22651 
19 2 1 3 2501 0.004 116.5 0.924 0.66 0.67 0.36 0.19 0.203 0.824 1.11960 
20 2 1 4 2480. 0.0039 111.6 0.932 0.63 0.70 0.37 0.18 0.212 0.806 1.10663 
21 2 2 1 2812.1 0.0050 119.2 0.822 0.80 0.66 0.32 0.24 0.199 0.773 1.13264 
22 2 2 2 2603 0.004 110.7 0.888 0.77 0.71 0.35 0.23 0.214 0.795 0.94908 
23 2 2 3 2503.0 0.004 104.8 0.923 0.76 0.75 0.36 0.22 0.226 0.845 0.83158 
24 2 2 4 2428.5 0.0043 100.4 0.951 0.71 0.78 0.38 0.21 0.236 0.799 0.80295 
25 2 3 1 3170.8 0.0049 108.3 0.729 0.79 0.73 0.29 0.23 0.219 0.705 1.24629 
26 2 3 2 2980 0.0041 100.6 0.775 0.66 0.78 0.31 0.19 0.236 0.709 1.27478 
27 2 3 3 2865.4 0.0036 95.29 0.806 0.59 0.83 0.32 0.17 0.249 0.597 1.25013 
28 2 3 4 2765.7 0.0031 91.24 0.835 0.50 0.86 0.33 0.15 0.260 0.581 1.27951 
29 2 4 1 3249 0.0042 102.3 0.711 0.68 0.77 0.28 0.20 0.232 0.558 1.41024 
30 2 4 2 3192.2 0.0038 92.19 0.724 0.62 0.85 0.28 0.18 0.257 0.578 1.34609 
31 2 4 3 3014.0 0.003 87.32 0.766 054 0.90 030 0.16 0.272 0.818 1.29628 
32 2 4 4 2809 0.0029 83.59 0.822 048 0.94 0.32 0.14 0.284 0.720 1.20629 
33 3 1 1 2894.4 0.0047 151.2 0.796 0.75 0.52 0.31 0.22 0.157 0.717 1.52240 
34 3 1 2 2735.2 0.0048 138.3 0.845 0.78 0.57 0.33 0.23 0.171 0.734 1.28541 
35 3 1 3 2641.9 0.0046 128 0.874 0.75 0.61 0.34 0.22 0.185 0.811 1.1832 
36 3 1 4 2520 0.0043 119.5 0.917 0.70 0.66 0.36 0.21 0.198 0.748 1.10086 
37 3 2 1 3305.1 0.0054 135.7 0.699 0.87 0.58 0.27 0.26 0.175 0.672 1.43900 
38 3 2 2 28132 0.0042 1243 0.821 0.68 0.63 0.32 0.20 0.191 0.733 1.39136 
39 3 2 3 2707.1 0.0041 115.1 0.853 0.67 0.68 0.34 0.20 0.206 0.727 1.25112 
40 3 2 4 2604.0 0.004 107.5 0.887 0.66 0.73 035 0.19 0.220 0.694 1.10704 
41 3 3 1 3514.1 0.0034 123.4 0.657 0.55 0.64 0.26 0.16 0.192 0.610 2.07097 
42 3 3 2 3390.7 0.0030 112.9 0.681 049 0.70 027 0.14 0.210 0.622 2.00608 
43 3 3 3 3313 0.0028 104.6 0.69 045 0.75 0.27 0.13 0.227 0.673 1.91873 
44 3 3 4 3137 0.0025 97.7 0.736 040 0.81 0.29 0.12 0.243 0.650 1.80251 
45 3 4 1 3721 0.0031 112.9 0.621 0.50 0.70 0.24 0.15 0.210 0.611 2.14538 
46 3 4 2 3481 0.0030 103.5 0.664 048 0.76 0.26 0.14 0.229 0.633 1.93399 
47 3 4 3 3249 0.0027 95.90 0.711 044 0.82 0.28 0.13 0.247 0.809 1.76709 
48 3 4 4 3025 0.0024 89.54 0.764 040 088 0.30 0.12 0.265 0.674 1.60353 
49 4 1 1 3150.5 0.0053 155.2 0.733 0.86 0.51 0.29 0.25 0.153 0.660 1.51310 
50 4 1 2 2961.5 0.0044 143.2 0.780 0.71 0.55 0.31 0.21 0.165 0.683 1.58673 
51 4 1 3 2785.7 0.004 135.5 0.829 0.67 0.58 0.33 0.20 0.175 0.756 1.49517 
32 4 1 4 2520 0.0040 129.4 0.917 0.65 0.61 0.36 0.19 0.183 0.740 1.27325 
53 4 2 1 3088 0.0050 139.5 0.748 0.80 0.56 0.29 0.24 0.17 0.67 1.48305 
54 4 2 2 2916 0.0042 128.9 0.792 0.67 0.61 0.31 0.20 0.18 0.68 1.51980 
55 4 2 3 2798.4 0.0040 121.9 0.826 065 0.64 0.33 0.19 0.194 0.694 1.42036 
56 4 2 4 2683.2 0.0036 116.4 0.861 0.58 0.68 0.34 0.17 0.204 0.669 1.3964 
57 4 3 1 3856.4 0.0033 127 0.599 0.53 0.62 0.23 0.15 0.187 0.574 2.3169 
58 4 3 2 3825.4 0.0030 117.2 0.604 049 0.67 0.24 0.14 0.202 0.583 2.2680 
59 4 3 3 3472.7 0.0027 110.7 0.665 0.44 0.71 0.26 0.13 0.214 0.639 2.1244 
60 4 3 4 3305.1 0.0026 105.8 0.699 0.42 0.74 0.27 0.12 0.224 0.618 1.9948 
61 4 4 1 3972.7 0.0030 116.2 0.581 0.49 0.68 0.23 0.14 0.204 0.583 2.3214 
62 4 4 2 3819.2 0.0027 107.3 0.605 044 0.73 0.24 0.13 0.221 0.590 2.2342 
63 4 4 3 3481 0.0026 101.5 0.664 0.43 0.78 0.26 0.12 0.234 0.510 2.0137 
64 4 4 4 3398.8 0.0024 97.00 0.680 0.38 0.81 0.27 0.11 0.24 0.272 1.9833 


Table 3 shows the effect of factors on measured re- 
sponse variables for weighting factor of w, = 0.4, wz = 
0.3 and w3= 0.3. The level which has the higher value 
when compared with other two levels is optimum levels 
for each factor. It is observed that the fourth level of fuel 
blends, injection pressure and first level of fuel tempera- 
ture with higher value when compared with other levels, 
and hence the levels (4, 4, 1) are taken as the optimum 
for assigned weighting factor of 0.4, 0.3, 0.3. 

Table 4 shows the effect of weighting factor on the op- 
timum combination levels of factors and percentage con- 
tribution of level influencing factors on the set objective. 
It is observed that the weighting factor plays an impor- 
tant role in deciding the contribution of factors on the set 
objective. As the weighting factor w; increases, the per- 
centage contribution of fuel blends on the set objective 
also increases which ensures the influence of fuel blends 
in spray tip penetration reduction. For the assigned w, the 
percentage contribution of fuel blends, injection pressure 
and fuel temperature depends upon the difference be- 
tween the weighting factors w and w3. As the difference 
between w. and w; for the same w, increases, the per- 
centage contribution of injection pressure and fuel blends 
increases while for fuel temperature it decreases till w, 
=0.4. If w, is more than 0.4, the increase in difference 
between w and w, for same w, results in increase in 
percentage contribution of fuel blends and injection 
pressure and decrease in percentage contribution of fuel 
temperature. 


ANOVA 


Table 4 shows the results of ANOVA for the weighting 
factor of w, = 0.4, w2 = 0.3 and w3 = 0.3. From the Table 
4, the percentage contribution of all the factors on the set 
objective can be observed. It is observed that fuel blends 
is the most influencing factor on the set objective since 
the percentage contribution is higher than the other two. 


Table 4 Results of ANOVA for w= 0.4, w = 0.3 and w3=0.3 


Fuel blends has a significant effect on the chosen objec- 
tive since its contribution is significant and the difference 
in percentage contribution between fuel injection pres- 
sure and fuel temperature also have some marginal effect 
as indicated in the ANOVA Table 5. Spray tip penetration 
depends on the fuel injection pressure and temperature. 
When compared with fuel blends, fuel injection pressure 
also has a considerable effect on all the response va- 
riables and its influence on the set objective is vital as 
obtained through ANOVA. It can also be observed that 
with this combination of weighting factors fuel tempera- 
ture has least effect on the objective since its percentage 
contribution is very low when compared with the other 
two. 

In the ANOVA analysis Table 5, there is a P value 
which is computed from F-ratio for each independent 
parameter in the model, if P value is less than 0.05, the 
parameters can be considered as statistically highly sig- 
nificant, since P- value for all the chosen parameters is 
less than 0.05, they are highly significant at 95% confi- 
dence level. From the analysis it is observed that fuel 
blends, injection pressure and fuel temperature are all the 
influencing factor for spray tip penetration, spray cone 
angle and SMD, since the value of P is zero in all cases. 
The purpose of analysis of ANOVA is to investigate the 
percentage contribution of variance over the response 
factor. Statistically, there is a tool called F-test named 
after Fischer to check the significance of variance on the 


Table 3 Effects of factor on response variables for w, = 
0.4,w2 = 0.3, w3= 0.3 


Factors Level 1 Level 2 Level 3 Level 4 
Fuel blends 0.869 1.017 1.59 1.80 
Injecti - 
A TAS 1.06 1.582 1.585 
sure (bar) 
Fuel tempera- 
1.45 1.40 1.32 1.257 


ture (°C) 


Percentage contribution of factors 


Combination Weighting factor Optimum level of factors from the ANOVA table 
No Fuel Injection Fuel Fuel Injection Fuel 
"i Wa di blends pressure (bar) temperature (°C) blends pressure (bar) temperature(?C) 

1 0.4 0.3 0.3 4 4 1 62.73 24.85 2.68 
2 0.4 0.4 0.2 4 4 1 48.94 40.01 0.11 
3 0.5 0.3 0.2 4 4 1 51.16 37.93 2.16 
4 0.5 0.4 0.1 4 4 2 39.065 51.68 0.098 
5 0.6 0.2 0.2 4 4 1 51.44 34.50 7.1 
6 0.6 0.25 0.15 4 4 1 46.16 42.76 3.84 
7 0.6 0.3 0.1 4 4 1 40.86 49.97 0.169 
8 0.7 0.2 0.1 4 4 1 41.18 47.02 4.4 
9 0.8 0.1 0.1 4 4 1 40.38 42 11.9 


Main Effects Plot (data means) for MRSN1 
Fuels 


Injection pressure 


Pe popa p 
oN R O o 
E T 7 


2 3 4 1 2 3 4 


Fuel temperature 


a 
œ 


Mean of MRSN1 
a 


1 2 3 4 


Fig. 2 Main effect plot for MRSN ratio 


output characteristics [13], usually when F>4 means the 
change of the design parameters has a significant effect 
over the output characteristics. 


Responses variables at optimized condition 


Confirmation experiment was conducted for the ob- 
tained optimum combination level and response variables 
of optimum combination have been compared with the 
response variables of the experiments conducted with the 
normal operating condition. (Fuel: KOME, Injection 
pressure: 240 bar, without preheating). Table 6 shows the 
response variables at optimized conditions and same is 
compared with the variables at normal operating condi- 
tion. It can be observed that spray cone angle is increased 
and spray tip penetration and SMD decreased as a result 
of this combined effect. The reason is attributed to im- 
proving spray pattern (due to its lower viscosity) and 
high rate of evaporation due to preheating. One of the 
optimum combination 4-4-1 is obtained by varying 
weighting factors and by Lg, full factorial design which 
shows an increase in spray cone angle and decrease in 


Table 5 Results of ANOVA for w\= 0.4, wz = 0.3 and w3= 0.3 


SMD and spray tip penetration. The fourth level of injec- 
tion pressure and fuel blends and first level of fuel tem- 
perature was the optimum combination for lower SMD 
and spray tip penetration and higher spray cone angle. 
This combination increases the spray cone angle by 33% 
with 2.7% and 0.8% decrease in SMD and spray tip pe- 
netration, when compared with the normal operating 
conditions. 


Conclusion 


This study explores the possibility of utilizing DOE 
and full factorial experimental design in spray characte- 
ristics of biodiesel and attempts to examine the combined 
effect of three factors on the three variables with the 
minimum number of experimental work. Experimental 
results were analyzed through statistical tools and find- 
ings of the analysis were used to make necessary deci- 
sion. In the present work, the optimum combination (4, 4, 
1) of fuel blends, injection pressure and fuel temperature 
in reducing SMD and spray tip penetration and increase 
in spray cone angle was obtained by MRSN. Different 
weighting factors were assigned to each response varia- 
ble to calculate the MRSN ratio and the obtained ratio 
was analyzed through ANOVA method. Fuel injection 
pressure and fuel blends appears to be most influencing 
factor in reducing SMD and spray tip penetration of the 
fuel next to fuel temperature with increase in spray cone 
angle. The conclusions are drawn with minimum number 
of experiments which saved time and cost because of the 
application of Taguchi based DOE. The optimum com- 
bination of factors was arrived after critical analysis of 
response variables with different possible combination of 
weighting factors. Without Taguchi method considerable 
experimental work has to carried out to study the effect 
of factors on the response variable. 


Factors DF Seq SS Adj SS Adj MS F P % Contribution 
Fuel blends(SS;) 3 8.6 8.66 2.88 745 0.0 62.73 
Injection pressure (SS+) 3 3.4 3.42 1.14 294 0.0 24.85 
Fuel temperature(SS,) 3 0.3 0.37 0.12 32.3 0.0 2.68 
Error (SS;) 27 0.1 0.10 0.00 
Total (SS;) 63 13 


Table 6 Effect of optimization on response variable 


Normal operating 


Optimized Condition 


% Variation of response variables 


Fuel Condition 
S e SMD Combination S e SMD S e SMD 
KOME 63.54 12.03 11.08 4-4-1 63.0101 18.01 10.78 0.8 33 2.7 
Diesel 55 12: 9.5 
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